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1. Introduction 

On. 

Feynman diagrams are the most important theoretical tool for particle 
physicists. They are an efficient link between theory and experiment. How- 
ever, their translation into actual numerical predictions is often very tedious 
if not impossible. Huge efforts have been devoted to their evaluation, and 
several powerful methods have been developed to systemize their treatment. 
The more complex a Feynman diagram is, the more important is it to find 
approximation procedures that allow to solve the problem with finite but 
reasonable accuracy. In this paper we will describe methods that have been 
developed over the recent years in order to systematically expand Feynman 
diagrams in their external parameters. 

2. Status of multi-loop calculations 

The complexity of a Feynman diagram with a certain number of loops 
mainly depends on its number of scales (i.e., masses and external momenta). 
In the one-loop case, the problem can be considered as solved. Any tensor 
integral can be reduced to integrals of unit numerator which have been 
studied extensively. Nowadays there are powerful software tools, on the 
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one hand concerned with the tensor reduction, on the other hand with the 
numerical or analytical evaluation of the integrals (see, e.g., [1]). 

At two-loop level the solution is not as general as in the one-loop case. 
However, the important class of two-point functions is well under control, 
and the development for three- and four-point functions is under continuous 
progress (see [2] for a list of references) . 

Therefore, two-loop calculations in theories like the electro-weak stan- 
dard model and even in supersymmetric models whose particle spectra give 
rise to Feynman diagrams with several different scales have become feasi- 
ble (e.g. [3]). 

Calculations at three-loop level mostly reside on two different classes of 
analytically solvable Feynman integrals. These two classes are 

• massless propagator-type diagrams where all internal lines are mass- 
less and only one external momentum is different from zero. Schemat- 
ically: 

t( \ f ad, a di a du P(q,h,k 2 ,k 3 ) 

I(n 1 ,...,n a ) = Jd fcxd k 2 d fc 3 (p 2 )wi .. ,(p2)n a ' (!) 

where a is the number of propagators and the pi are linear combina- 
tions of the kj and the external momentum q. 

• massive tadpole diagrams not carrying any external momenta and 
internal lines being either massless or carrying a common mass m. 
Schematically: 



j(m,... ,n„)= / d D hd D k 2 d D k 3 / 2 . :;„: i\ ^ ■ ( 2 ) 



,i>;. Alh. aL>u P(h,k 2 ,h) 

(ml + pj) ni ■ ■ ■ (m 2 a + p 2 a ) n ° 

where the rrii are either equal to zero or m, and the pi are linear 
combinations of the loop momenta kj . 

P(. . . ) is a polynomial of products of its arguments. The method how to 
solve such integrals is called the integration- by-parts algorithm [4]. It is 
based on identities derived from the fact that the D-dimensional integral 
over a total derivative is equal to zero: 

d D p^-f(p,...)=0. (3) 

These identities can be arranged in such a way that they yield recurrence 
relations that allow to reduce some of the "indices" m, ... , n a in (1), (2) to 
zero. At three-loop level, these relations have been derived in [4] for massless 



propagators and in [5,6] for massive tadpole integrals. Their application to 
a general three-loop diagram may generate huge intermediate expressions 
that easily exceed several hundreds of megabytes on a computer. This is 
why one needs to implement the relations to powerful computer algebra 
systems like FORM or REDUCE. Two such implementations are MINCER [7], 
concerned with the massless propagator diagrams and MATAD [8], dealing 
with the massive tadpoles (see [2] for a review on automatic computation 
of Feynman diagrams) . 

The two classes of single-scale diagrams mentioned above already have 
a huge number of important applications. The most popular one probably 
is the total cross section for hadron production in e + e~ annihilation (see, 
e.g., [9]), usually written as the hadronic R ratio, in the limit of vanishing 
quark masses. Using the optical theorem, it can be expressed through the 
imaginary part of the photon polarization function: 



R(s) = 127rlml% 



2- 



(4a) 



q 2 =s+ie 






where I% 2 ) = \ 1?"™* U^q) , (4b) 



and IV(g) = * / d 4 xe^ x (0\Tj^x)j u (0)\0) , ^ = jry^ ■ ( 4c ) 



ip is a quark field of mass m. The diagrams contributing to IL^ (q) up to two- 
loop order are shown in Fig. 1. Up to three-loop order, such diagrams can 
directly be computed by the program MINCER mentioned before {nota bene 
in the massless limit!). Another important application is the computation 
of moments of the polarization function, d n /d(q 2 ) n ■ H(q 2 )\ q 2 = Q. They can 
be obtained by applying the derivatives and the nullification of q 2 before 
performing the loop integrations [10]. Furthermore, since renormalization 
group functions like the QCD (3 function or anomalous dimensions in the 
MS scheme are independent of any masses and momenta, their evaluation 
can be performed by computing single-scale diagrams. 

However, returning to the R ratio defined above, the limit of vanishing 
quark mass may not be satisfactory, especially if one is interested in energy 
regions not too far above one of the quark thresholds. As long as the exact 
evaluation of three-loop diagrams involving a non- vanishing mass as well as 
an arbitrary external momentum is not possible, one may hope to reduce 
the integrals to single-scale diagrams by performing an expansion in the 
quark mass. In the optimal case, a finite number of terms in the expansion 
will approximate the full result to reasonable accuracy, and the inclusion 
of higher order terms will gradually decrease the error. To get an idea on 
what the result should look like, let us consider the exact one-loop result 
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Fig. 1. Diagrams contributing to the photon polarization function 11(g). The outer 
line denotes a photon of momentum q, the plain lines are quarks and the spiral 
ones denote gluons. 



for the photon polarization function: 
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and jit is the renormalization scale. The pole 1 in e = {D — 4)/2, where D is 
the space-time dimension, will eventually disappear upon global renormal- 
ization (e.g., by requiring n(0) = 0). The result of (5) can be expanded in 
terms of small mass m, yielding: 
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(6) 
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with l qili = ln(—q 2 /fj, 2 ) and l qm = ln(—q 2 /m 2 ). One observes that the co- 
efficients of the series in m 2 jq 2 contain non- analytical pieces in terms of 
logarithms. They develop an imaginary part by means of 



ln(— s — ie) = —m + In s for s > 0. 
The hadronic R ratio is therefore given by (see (4a)): 

,2\2 \ 



R(s) 
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+ ... 



(7) 



(8) 



The accompanying ln47r and 7e are suppressed throughout the paper. 



The question is now if it is possible to obtain the expansion given in (6) 
directly from the Feynman integrals, i.e. without having to know the exact 
result. As a first guess one may try to perform a Taylor expansion of the 
integrand, thereby arriving at massless propagator diagrams. However, it 
is clear that such a "naive Taylor expansion" can not be the whole answer. 
For example, it is impossible to reproduce the logarithmic mass dependence 
in this way. Nevertheless, let us look at the result: 
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In fact, the first two orders in m 2 /q 2 are reproduced correctly. The m 4 /g 4 
term, however, is completely different, and there is even an additional pole 
in e. Only the logarithmic q 2 dependence is reproduced. Thus, in general 
the naive Taylor expansion is not sufficient to arrive at the desired result. 
However, in the next section we will see that by including well-defined ad- 
ditional terms one indeed can obtain the correct expansion. 

3. Asymptotic Behavior 

This section is divided into three parts, all concerned with the problem 
of expanding Feynman diagrams in their external parameters, as it was 
raised in Section 2. Cross-references between the three parts of this section 
will demonstrate the close correspondence of the individual formulations. 

In Section 3.1, the problem will be approached from a field theoretical 
point of view. The resulting expansion will be derived from the operator 
product expansion formulated in the MS scheme. The viewpoint of Sec- 
tion 3.2, on the other hand, examines the individual Feynman integrals 
that contribute to a certain problem. By a thorough investigation of the in- 
tegration regions for the loop momenta and a subsequent Taylor expansion 
in the appropriate variables, one can derive rules that allow to obtain the 
expansion of the full result in a very efficient way. 

In certain cases these rules could be phrased in a mainly diagrammatical 
language. For diagrams involving large external momenta or large masses 
this graphical formulation will be described in Section 3.3. 

3.1. Operator Product Expansion 

The asymptotic behavior of the two-point correlator of (4c) in the limit 
—q 2 = Q 2 — ► oo is formally known to all orders of perturbation theory. It 
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Fig. 2. Sample diagrams contributing to (Oi) and (C2) 



is given by an operator product expansion (OPE): 

U^(q) =iJd 4 xe^ x (0\Tj ll (x)j u (0)\0) " 9 ^°° £C n)/iI/ (0 n ) . (10) 



The C n ^ u are complex functions, and the O n are operators composite of 
fields of the QCD Lagrangian. We keep only Lorentz scalar operators 
because all others vanish when sandwiched between the vacuum states. 
Transversal coefficient functions C n will be defined in analogy to the Eq. (4b). 
The operators are usually sorted according to their mass dimension. It is 
convenient to allow only operators of even mass dimension which is achieved 
by appropriate factors of the quark mass m (only one quark shall be consid- 
ered as massive for the sake of clarity). Up to dimension four, the following 
set of operators is relevant: 

0(°) = 1 , 0( 2 ) = m 2 , (11) 

where G^ u is the gluonic field strength tensor and ip is again the quark 
field (the superscript "(4)" of the dimension-4 operators will be dropped in 
what follows). If the operators are understood to be normal ordered, the 
vacuum expectation values of the non-trivial (i.e. not proportional to unity) 
operators are equal to zero in perturbation theory. In such an approach these 
operators are used to parameterize non-perturbative effects (see, e.g., [11]). 

On the other hand, if one abandons normal ordering and applies mini- 
mal subtraction, the vacuum expectation values of the field operators receive 
also perturbative contributions. The corresponding diagrams are massive 
tadpoles which by definition only depend on the quark mass and the renor- 
malization scale [i. Examples are shown in Fig. 2. They lead to the following 



results: 
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(12) 
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Only the first non- vanishing order in a s is quoted here. However, the results 
for the operators that are proportional to unity (m 2n ,n = 0, 1,2) are valid 
to all orders of perturbation theory by definition. Note that there are still 
poles in e which is why the operators are marked with the superscript B. 
These poles disappear upon global renormalization, thereby inducing mixing 
of the operators according to 

O n = J2 Z nmO*. (13) 



For the dimension-4 operators the renormalization matrix Z nm in the MS 
scheme was computed in [12] by expressing it in terms of the charge and mass 
renormalization constants of QCD, plus the one for the QCD vacuum energy. 
They are currently known to 0(aj.) [13,14] and O(o%) [15], respectively. The 
results for the renormalized vacuum expectation values of 0\ and 02 U P to 
(D(a s ) can be found in [16]. 

The crucial observation about using the MS scheme was made in [17]: 
It appears that in this approach the coefficient functions are independent 
of the quark masses. They only depend on the external momentum q and 
the renormalization scale [i. It was shown that their computation can be 
reduced to the evaluation of massless propagator- type diagrams. The most 
refined method for this purpose is called the "method of projectors" [18]: 
let us define "bare" coefficient functions through 

Y,C n , B 0* = Y,CnO n . (14) 

n n 

Then the ones up to dimension four, for example, are obtained in the fol- 



lowing way (only sample diagrams are displayed here) [19]: 
[C B 0) , C B 2) , C 3iB ] = 



d 1 d 2 




me=0 



m,B=p=0 

(15) 

where p is the momentum carried by the external quark-, gluon-, and ghost- 
lines (the latter arise only at higher orders in a s ). V\ and V2 are "projectors" 
depending polynomially on the derivatives w.r.t. p and me- For example, 



V 2 



d d 



dp 1 dms 



1 m ( d 1 v d 

Tr 1 Y 

4n c \drriB D dp h 



(16) 



It is understood that the derivatives act on the integrands and the nulli- 
fication of p and m-Q is performed before integration. Note, however, that 
the momentum carried by the external currents (wavy lines), is q 7^ 0. Ob- 
viously, C B , C B , and 63,3 are just the coefficients of m B , m B , and to b 
of the naive Taylor expansion, respectively. The expressions for Ci B and 
C2 Bj on the other hand, read as follows: 
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where we have kept the terms up to 0(e) because they contribute to the 
finite part of II(g 2 ). In fact, with these ingredients it is possible to compute 
the polarization function up to order a®m 4 . Diagrammatically one finds 
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(18) 



where in the second term on the right hand side the projection with Vi 
from (16) is implicit. The factor 2 arises from the symmetrical diagram 
that also contributes to C2,b- Note that there is no contribution from 0\ 
at this order. The Taylor expansion in the first term is to be carried out up 
to njg. This first term is given by Eq. (9), since m = ttt-b + 0{a s ). Using 
the results of (12) and (17) for C^b and (Of), one obtains the result for 
the second term of (18): 
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(19) 



Adding it to (9) exactly reproduces the result for H(q 2 ) up to 0(m 4 ) given 
in (6). 

The conclusion from these considerations is that in addition to the 
naive Taylor expansion of (9) one should include an extra term, given by 
2C2 i b(02 3 )' i n or der to arrive at the correct result for the polarization func- 
tion. However, the important point about Eq. (18) is that the original 
diagram is reduced to single-scale factors (i.e., massless propagators and 
massive tadpoles). 



3.2. Strategy of Regions [20,21] 

Let us for the moment forget about OPE again and consider only the 
Feynman integral for the one-loop diagram, concentrating on the scalar case 
for the sake of clarity: 
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m 2 + k 2 m 2 + [k — Q) 2 ' 



(20) 



where the momenta are taken in Euclidean space and integration is over k. 
Assume now that m 2 <C Q 2 (= — q 2 ). The integral may be split into the 
following regions: 

(i) : k 2 > m 2 and (k - Q) 2 > m 2 

(m) : k 2 ~m 2 => (k- Q) 2 » m 2 (21) 

(iii) : (k - Q) 2 ~ m 2 => k 2 » m 2 , 



where ~ means "of the order of". In region (i), the integrand can be 
expanded in terms of small m: 

1 1 f _ 1 1 



('■) 



m 2 + k 2 m 2 + (k — Q) 2 



J. jr 



('■) 



k 2 (k - Qf 



m 
1? 



(0 



+ ... 



m 2 + k 2 m 2 + (k — Q) 2 

2 x (22) 
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In region (ii), k is considered to be of the same order of magnitude as m, 
so one should expand in m and k at the same time: 



/, 



i i r i 

1 m,k 



m 2 + k 2 m 2 + {k — Q) 2 Ju^ ' m 2 + k 2 m 2 + (k — Q) 2 
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1 1 ( m 2 + k 2 -2k-Q 

J- "T7J l - • • • 



m 2 + k 2 Q 2 V <5 2 



Region (Hi) can be mapped onto region (ii) by substituting k' = Q — k, 
meaning k' 2 ~ m 2 =4* (k 1 — Q) 2 3> w 2 , and we arrive at an expression 
analoguous to (23). 

Now we add and subtract the complementary regions to the integrals 
above and find: 

1 f 



m 2 + k 2 m 2 + (k — Q) 2 



r _J I + 2 f 1 t k I ( 24 ) 



'm 2 + A; 2 m 2 + (/c — Q) 2 j m 2 + /c 2 ' m 2 + (k — Q) c ' 

C, 



where 



r / t 1 x i r ' 



/ 



(ii)u(m) ™ 2 + k 2 m 2 + {k- Q) 2 J^ui) m 2 + /c 2 
1 /" 1 _ 1 

(i)U(ii) "Z 2 



(25) 



X T »^ m 2 + {k _ Q) 2 + A m 2 + fcA* m 2 + {k , _ Q y 



In each of the different regions, one can again expand w.r.t. the appropriate 
parameters, for example: 



(jj) m 2 + /c 2 m 2 + (k — Q) 2 

{ii) K Tm m 2 + k 2) [ T ^ m 2 + {k _QY 



(26) 



Finally, for C one finds: 



i(i)u(ii)u(m) V m "i 2 + fc 2 y V m ' k m 2 + (k - Q) 2 J 

which corresponds to massless tadpole integrals and therefore vanishes in 
dimensional regularization which we are using throughout. 2 

2 I acknowledge a useful conversation with K. Melnikov on this issue. 
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The second term on the r.h.s. of Eq. (24) can be represented graphically 



as 



2 • ^ W, (28) 

meaning that the solid line should be expanded in all quantities but Q. 
Doing this, the actual integral will be of the tadpole type, with a more or 
less complicated vertex insertion. It is obtained from the original diagram 
by shrinking the solid lines shown in (28) to a point. In this sense, one can 
again use the diagrammatical equation (18) in order to represent Eq. (24). 
However, the interpretation is slightly different: In the second term, the 
diagram left of "*" means the solid line of (28), and the one right of "*" is 
the tadpole which remains when shrinking the expanded lines to a point. 
But the main difference between both interpretations of Eq. (18) is that 
according to Section 3.1 it is only valid up to terms of order (m 2 /q 2 ) 2 ; 
according to Section 3.2, however, one may obtain the diagram to the r.h.s. 
up to arbitrary high powers in m 2 /q 2 by pushing the Taylor expansions on 
the l.h.s. to sufficiently large order. 

3.3. Large momentum and hard mass procedure [17, 22, 18, 23-25] 

So far, two different viewpoints for the small-mass expansion of the 
polarization function have been presented. One is based on OPE and uses 
the language of quantum fields, while the second one examines the very 
Feynman integrals. The third viewpoint we will focus on is formulated 
mainly in terms of Feynman diagrams. Given a certain Feynman diagram 
with a particular distribution of masses and external momenta, the method 
generates a set of simpler diagrams which correspond to the asymptotic 
form of the original diagram in certain limits of the external parameters. 

The particular case of a large external momentum as it was considered 
so far is usually referred to as "large momentum procedure" . The opposite 
case of a mass being much larger than any other scale is called "hard mass 
procedure" and will be addressed below. 

The prescription for the asymptotic expansion of Feynman diagrams can 
be summarized by the following formula [25]: 

F{T) - X>(r\ 7 ) *T^( 7 ). (29) 

7 

Here, T is the Feynman diagram under consideration, and J-(T) is the cor- 
responding Feynman integral. It shall contain either a set of large external 
momenta {Q} or of large masses {M}. The arrow (— ►) denotes that the 
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r.h.s. is valid in the asymptotic limit of the M or Q going to infinity. The 
sum goes over all subgraphs 7 of T that fulfill certain conditions to be de- 
scribed below. r\7 means the diagram that results when, within T, all 
lines of 7 are shrunk to points. T means Taylor expansion w.r.t. all masses 
and external momenta that are not large. In particular, also those exter- 
nal momenta of 7 that appear to be integration momenta in T have to be 
considered as small. The Taylor expansions are understood to be applied 
before any loop integrations are performed. In the following we will refer to 
the 7 as hard subgraphs or simply subgraphs, to T\7 as the corresponding 
co-subgraphs. The "•" means that T^F(j) shall be inserted into J r (T\7) at 
the point to which 7 was contracted. 

In other words: within I, all propagators of 7 have to be expanded w.r.t. 
the masses and external momenta of 7 that are not large. 

3.3.1. Large momentum procedure 

For the particular case of the large momentum procedure, the conditions 
that specify the hard subgraphs are as follows: 

Every 7 has to (i) contain all vertices where a large momentum enters 
or leaves the graph and (ii) be one-particle irreducible if these vertices were 
connected by an extra line. 

As an example, consider again the one-loop diagram contributing to 
the photon polarization function, Fig. 1(a), in the limit of large external 
momentum. The set of hard subgraphs that emerge consists of three dia- 
grams: first there is the diagram itself; the corresponding co-subgraph is 
just a point. The second subgraph is the one shown in (28) if the dashed 
line is omitted, and the third subgraph is the one symmetrical to that. The 
corresponding co-subgraphs are one-loop tadpole diagrams. In this way one 
again arrives at Eq. (18). The interpretation of the terms is the same as it 
was in the approach of Section 3.2. 

As a two- loop example, let us examine the diagram shown in Fig. 1 (b): 
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(30) 



13 



where the subgraphs are given by the solid lines of the diagrams left of "*" . 
The first term on the r.h.s. corresponds to the naive Taylor expansion. Note 
that the last diagram is zero in dimensional regularization. 

It is interesting to relate this set of terms to the viewpoints of the pre- 
vious section. In the operator language, the first term corresponds to the 
trivial operators m 2n . The second one (actually only its m 4 /q 4 contribu- 
tion) is related to Cip,{0^ ), the last one to C\^{Of). The other two terms, 
however, have no correspondence to any of the operators of Section 3.1. This 
is due to the fact that we considered only operators up to dimension four. 
But relation (30) is valid up to arbitrary orders in m 2 /q 2 . So the conclu- 
sion is that the third and fourth term on the l.h.s. are of order (m 2 /q 2 ) 3 or 
higher. 

The viewpoint described in Section 3.2, on the other hand, reproduces 
exactly the same terms as shown above. 



3.3.2. Hard mass procedure 

So far only the case of an external momentum being much larger than 
any other scale of the problem was considered. In this section the so-called 
hard mass procedure will be discussed. Here it is assumed that the diagram 
carries a mass that is much larger than all other masses and external mo- 
menta. Equation (29) remains valid, only the classification of the subgraphs 
is different. In the case of the hard mass procedure, 7 must (i) contain all 
lines carrying a large mass (ii) be one-particle irreducible in its connected 
parts after contracting the heavy lines. Consider the following two-loop 
diagram as an example: 



q+k+l^c^ ""Nj q+k 

(31) 



ijx 



where q is the external momentum flowing through the diagram from right 
to left. The mass of the thick line will be denoted by M. The imaginary part 
of (31) contributes to the electro- weak one- loop corrections of the process 
Z — > bb if top-quarks are attributed to the thick lines, b quarks to the plain 
thin lines, W bosons to the inner, and Z bosons to the outer wavy lines (see 
also Section 4 below). In the limit q 2 <C M 2 , the following subdiagrams 
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emerge: 



k 2 ~M 2 
I 2 ~M 2 

k 2 ~M 2 
I 2 <M 2 

k 2 «M 2 



• vA^V^ , 



(32) 



I 2 



M 2 



•k -a* W j 



fc 2 «M 2 
Z 2 <M 2 



*■- * 




where we have indicated the region of loop momenta that generates the 
corresponding subdiagram according to the considerations or Section 3.2. 
It is again understood that the solid lines are to be expanded w.r.t. all 
external momenta and masses except M. The co-subgraphs, shown right of 
"*", are obtained by contracting the solid lines to points. 

3.3.3. Successive application of large momentum and hard mass procedure 

In theories with many different particles of various masses, a realistic 
process generally involves several scales. The successive application of the 
large momentum and the hard mass procedure can be used to reduce any 
Feynman diagram to single-scale factors in this case, as long as a hierarchy 
among the different scales can be defined. This strategy was followed, for 
example, in the calculation of the anomalous magnetic moment of the muon 
[26], and later, using an automated setup, to the decay of the Z boson into 
b quarks [27] (see Section 4 below). 



3.4- Expansions in other limits and range of applicability 

Two concluding remarks shall be made to complete this formal section. 
First it should be noted that asymptotic expansions have been developed for 
various limiting cases (see, e.g., [28,20]). Only the two most straightforward 
ones have been described above. Other situations arise, for example, if the 
external momentum is either close to a mass or a threshold of the diagram. 
It appears that the most convenient way to formulate these expansions is in 
terms of the language of Section 3.2. Applications will be discussed briefly 
in the next Section. 

The second remark that should be made concerns the applicability of 
asymptotic expansions, or in other words, the convergence properties of 
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the series. By investigating the global structure of the problem, one often 
can read off regions of convergence for the function to be approximated. 
These regions may extend to values of the expansion parameter that are 
way beyond the initially required conditions. 

Restrictions on the region of convergence, on the other hand, are mostly 
induced by the presence of thresholds. Thus, it often appears that with 
increasing number of loops the range of validity for the result decreases 
due to additional thresholds that are absent at lower order of perturbation 
theory. 

4. Asymptotic expansions in practical applications 

R(s) to 0(a 2 i ): in [29] the method of Section 3.3.1 was directly applied 
to the three-loop polarization function. For this purpose, the diagrammati- 
cal prescriptions were implemented in a computer program, called Imp [30], 
in order to cope with the large number of subdiagrams that had to be gen- 
erated. For details on the calculation and the discussion of the results and 
their convergence properties we refer to [29,31]. 

RG functions in MOM scheme: the large momentum procedure was 
also used to evaluate the relations between the renormalization constants 
of QCD in the MS scheme to the ones in other schemes, in particular the 
MOM scheme, at 0(a^). See [32] for more details. 

Z — > b6 to order Gfcz s : a sample diagram whose imaginary part 
contributes at 0{G-p) is shown in Eq. (31). At order Gf« s , an additional 
gluon has to be attached. One can apply the hard mass procedure w.r.t. 
Mi in this case, and it turns out that the leading term is proportional 
to M 2 due to the large mass splitting between the bottom and the top 
quark. There are two other scales in the problem, M\\r and Mz- In [27] 
they were factorized by a successive application of the hard mass procedure 
according to My/ 3> Mz-, an inequality whose use clearly has to be justified 
(see [27,33]). The whole calculation would not have been possible without 
the computer program EXP [34]. As compared to Imp, it also performs 
the hard mass procedure and its combination with the large momentum 
procedure. Meanwhile the results of [27] have been confirmed by a different 
approach [35]. 

t — >■ bW to O(a^): the problem here is that the contributing three- 
loop diagrams are actually on-shell, q 2 = M 2 . Two different methods have 
been used to compute this process: The first calculation [36] used asymptotic 
expansions in the limit 1 — M 2 jM 2 <C 1, the second one [37] evaluated the 
off-shell diagrams for q 2 <C M 2 with the help of the hard mass procedure, 
taking the limit q 2 — > M 2 after performing a Pade approximation (see [37, 
38] for details). The perfect agreement and the fairly high accuracy of 
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the results in both approaches is a clear demonstration of the power of 
asymptotic expansions. 

e+e — — ► tt near threshold: the development of asymptotic expan- 
sions in the threshold limit allowed to compute this process up to 0(v 2 , 
a s v, a 2 ) [39], where v is the velocity of the top quarks in the cms system. 
It turned out that the corrections are huge and exhibit a strong depen- 
dence on the renormalization scale. The interpretation of this result and its 
implications are still an ongoing discussion (see [40] and references therein). 

MS to pole mass conversion at O(a^): among other things, the 
progress in threshold calculations mentioned above makes it very important 
to be able to express the MS mass in terms of the pole mass at O(a^). The 
corresponding calculation of this relation was based on asymptotic expan- 
sions in the limit of small and large quark mass and interpolated the result 
for the actual on-shell diagrams with the help of Pade approximations [41]. 

/i. decay and semileptonic b decays: for special final state config- 
urations the decay b — > civ was computed up to 0(a 2 ) using asymptotic 
expansions in the limit 1 — M 2 /M 2 <C 1 [42]. For b — ► ulv, the full inclu- 
sive result was obtained through a four-loop calculation, where the actual 
integrals were calculated by performing an asymptotic expansion and re- 
summing the full series [43]. The same technique was applied to the 2-loop 
QED corrections to muon decay [44]. Using an approach similar to the one 
of [37] (see above), the results of [44,43] were recently confirmed [45]. 

5. Back to OPE: RG improvement 

The advantages of the approach of Sections 3.2 and 3.3 in comparison 
to the one of Section 3.1 when calculating Feynman diagrams are obvious: 
In the OPE language one needs to define a full set of operators for every 
order in m 2 and to apply the appropriate projectors in order to obtain the 
coefficient functions. Using the large momentum procedure, all operations 
can be performed on the level of Feynman diagrams. The depth of the 
expansion in m? is just a matter of evaluating the Taylor expansions to 
sufficiently high orders. However, the loss of contact to the quantum field 
theoretic level has its price. In this section two advantages of the OPE 
language over the strict application of the large momentum procedure will 
be described. 

5.1. Resummation of logarithms [46] 

Consider the expansion of the polarization function w.r.t. small m 2 /q 2 
(see Eq. (6)). The coefficients of this series contain logarithms of the form 
ln(^ 2 /m 2 ) and ln(— fi 2 /q 2 ). The former arise from the tadpole diagrams (in 
the OPE language: the vacuum expectation values; in the large momentum 
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procedure: the co-subgraphs), the latter from the massless propagator di- 
agrams (coefficient functions/hard subgraphs). Working at fixed order in 
perturbation theory, one has to choose some value for [i in order to make 
predictions for physical quantities. It is clear that one should choose fj, such 
that the coefficients of the perturbative series are small. However, with the 
two types of logarithms above this might be impossible. In the OPE ap- 
proach this problem can be solved in the following way: the operators fulfill 
the renormalization group (RG) equation 

" n 

where j nm is their anomalous dimension related to the renormalization ma- 
trix Z nm of (13) by 

Inm = 2_^ ( M -r^2 Z mk J (Z~ )kn • (34) 

The solution of (33) is of the form 

On(fJ-) = ^ R nmO m {no) , (35) 



where R nm depends on \x and \xq only implicitly through its dependence on 
a s . One can now rewrite the OPE as 

I% 2 ) = Y^ C n (jj)O n (ji) = Y C n (fi)R nm O m (ji ) . (36) 



Setting fj? = q 2 and /j,q = m 2 removes both types of logarithms discussed 
above. Using the diagram-wise expansions of Sections 3.2 or 3.3, such kind 
of resummation is not possible in a straightforward way. 

5.2. Reconstruction of higher orders in a s 

Consider the expansion of Ti{q 2 ) in small m 2 /q 2 at (I + l)-loop level. It 
appears that the only proper (Z+l)-loop diagrams that contribute arise from 
the naive Taylor expansion, or in other words, from the trivial operators 
m 2n . All the additional terms are products of diagrams with a lower number 
of loops. In the following we will show that this fact allows to derive the 
physically relevant quantity R = Imll at (I + l)-loop level from an £-loop 
calculation assuming the OPE approach of Section 3.1. 
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To be specific, consider the calculation of the m 4 terms of U(q 2 ). The 
dimension-4 piece of the OPE fulfills the following RG equation: 

i 3 
n=\ 

The effect of the derivative on the O n is again determined by (33) . On the 
other hand, because the C n do not depend on masses, we can write 

f£?°'-(ik +aJ> &) c - (38) 

where L = ln(—fi 2 /q 2 ). Inserting this into (37), collecting the coefficients 
of 03 and using 733 = 47™" (this can be seen by recalling fi 2 -^m = j m m), 
one obtains [16]: 

— C 3 = -4 7 m C 3 - a s f3—C 3 - Y, 7n 3 C„ . (39) 

° s n=l,2 

Performing an /-loop calculation, C\ and Ci are known to order a s , C3 only 
to a^ -1 . But since (3 and 7™ do not contain terms of order 0^, the r.h.s. 
is known to 0{a s ) (assuming that the anomalous dimensions are known to 
appropriate order), and so is the l.h.s. The logarithmic terms (and thus the 
desired imaginary part) of C3 can therefore be obtained to 0(a s ) by trivial 
integration. 

This strategy was first followed in [16] to derive the m 4 a 2 terms of R(s). 
The extension to m^a^ was performed in [47] (see also [48]). 

6. Conclusions 

In this lecture we have discussed the methods for asymptotic expansions 
of Feynman integrals. The field theoretical approach via operator product 
expansion has been compared to the diagram-wise methods, and the advan- 
tages of both strategies have been outlined. We hope that the importance 
of the numerous applications — we could sketch only a few of them — has 
convinced the reader of the power and flexibility of these expansions. 
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